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Abstract 

QCD at finite isospin chemical potential /i/ is studied. This theory has no fermion 
sign problem and can be simulated on the lattice using present-day techniques. We 
solve this theory analytically in two limits: low /ij where chiral perturbation theory 
is applicable, and asymptotically high /// where perturbative QCD is at work. At 
low isospin density the ground state is a superfluid pion condensate. At very high 
density it is a Fermi liquid with Cooper pairing. The pairs carry the same quantum 
numbers as the pions. Motivated by this observation, we put forward a conjecture 
that the transition from hadron to quark matter is smooth. The conjecture passes 
several nontrivial tests. Our results imply a nontrivial phase diagram in the space of 
temperature and chemical potentials of isospin and baryon number. At asymtotically 
large values of m and small values of baryon chemical potential the ground state is 
in a phase similar to Fulde-Ferrell-Larkin-Ovchinnikov phase. It is characterized by 
a spatially modulated superfluid order parameter (uj^d) and may be the asymptotic 
limit of the inhomogeneous pion condensation phase advocated by Migdal and others. 



1 Introduction 

Good knowledge of QCD in the regime of finite temperature and baryon density is crucial 
for understanding a wide range of physical phenomena. In cosmology, one faces the problem 
of understanding how the Universe has evolved through the QCD phase transition at tem- 
perature T ~ 150 MeV. Due to the smallness of the baryon asymmetry, finite-temperature 
QCD should be sufficient to deal with this problem. However, for the physics of heavy-ion 
collisions, one needs to know how QCD behaves when both temperature and baryon chemical 
potential are finite. Lastly, neutron stars require the knowledge of matter in the "dense" 
regime, i.e. at large baryon densities and very low temperatures. Much less is known about 
the last two regimes compared to that of high-temperature baryon-antibaryon-symmetric 
QCD. 
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Let us use neutron stars as an example to illustrate the range of questions one would like to 
have answers to. The equation of state (EOS) of nuclear matter at high densities determines 
the mass-radius relationship and the maximum mass of neutron stars. Walecka model of 
nuclear matter predicts that the EOS becomes stiffer at higher densities and approaches the 
Zel'dovich limit, e = p (where the velocity of sound approaches light speed) at very high 
densities. However, at asymptotically high densities one expects nuclear matter to become 
a weakly-interacting quark liquid, with a much softer EOS e = 3p. At what density does the 
transition happen, and is it a phase transition or a crossover? 

Migdal and others 0] suggested that at very high densities pion condensation might 
happen. It is also argued that at even higher densities kaons are condensed 0. One would 
like to know whether pion and kaon condensations do indeed occur in nuclear matter, before 
the transition to quark matter has happened. 

Finally, there is a strong recent interest in the phenomenon of color superconductivity 
0, ^. One very interesting prediction is that at high enough chemical potential, the ground 
state of QCD is the "color-fiavor-locking" state @, which breaks chiral symmetry. However, 
while reliable results can be obtained at asymptotically high densities, where the strong 
coupling is small 0, § |^, it is not known how to extend these results to the region of 
smaller, more realistic, densities without having to rely on uncontrollable approximations. 

Lacking reliable analytical means to approach QCD in the strong coupling regime, one 
naturally turns to numerical methods. First principle lattice numerical Monte Carlo calcula- 
tions provide a solid basis for our knowledge of the finite temperature regime. However, the 
regime of finite baryon chemical potential fiB is still unaccessible by Monte Carlo, because 
present methods of evaluating QCD partition function require taking a path integral with a 
measure which includes a com])/ea; fermion determinant. At zero chemical potential one can 
simply ignore the determinant (as one does in the popular quenched approximation) and still 
find reasonable results for physical quantities. However, at finite fiB this procedure leads to 
qualitatively unacceptable answers as was realized long time ago |T^. It has been under- 
stood more recently that the quenched approximation breaks down at finite /i^ because it 
describes an unphysical theory containing, beside the normal quarks, the so-called conjugate 



quarks with opposite baryon charge [11 



As a side remark, one notes that though the conjugate quarks are absent in real QCD, 
there are many theories where they are naturally present. One class of such theories is QCD 



with two colors, where quarks are self-conjugate [12|. Another class contains theories with 



quarks in the adjoint color representation |T^. In all these theories, the positivity of the 
fermion determinant ensures the applicability of lattice Monte Carlo methods. However, the 
particle content of all these theories is very different from the real world. 

The failure of the quenched approximation in real QCD at finite fis and our inability to 
include complex fermion determinant in a Monte Carlo simulation is one of the main reasons 
for our understanding of QCD at finite baryon density to be still rudimentary. 

One way QCD at finite baryon density is different from finite-temperature QCD is that 
the transition from hadronic to quark degrees of freedom occurs due to the large density of 
a conserved charge (such as baryon number) while temperature plays no role. This is the 
motivation for us to turn to QCD at finite chemical potential /i/ of isospin (more precisely, 
of the third component of isospin, J3), which is conserved by strong interaction. 



2 



Before going into details, we would like to comment on the relevance of this regime to 
the real world. Nature does provide us with non-zero fij systems in the form of isospin- 
asymmetric matter (e.g., inside neutron stars), however, the latter contains both isospin 
density and baryon number density. In contrast, the idealized system considered in this 
paper does not carry baryon number: the chemical potentials of the two light quarks, u and 
d, are equal in magnitude, 1/1/1/2, and opposite in sign. Such a system, strictly speaking, is 
unstable with respect to weak decays which do not conserve isospin, and, as we shall see, is 
also not electrically neutral and thus does not exist in the thermodynamic limit. However, 
since we are interested in the dynamics of strong interaction alone, one can imagine that all 
relatively unimportant electromagnetic and weak effects are turned off. Once this is done, 
we have a nontrivial regime which, as we shall see, is accessible by present lattice Monte 
Carlo methods, while being analytically tractable in various interesting limits. As a result, 
the system we consider has a potential to improve substantially our understanding of cold 
dense QCD. This regime carries many attractive traits of two-color QCD |]12|, |13|, but is 
realized in a physically relevant theory — QCD with three colors. 



2 Positivity and QCD inequalities 

Since in Euclidean space the fermion determinant of our theory is real and positive, some 



rigorous results on the low-energy behavior can be obtained from QCD inequalities p^ , |T3 
Let us recall how the inequalities are derived in vacuum QCD. The starting point is the 
following property of the Euclidean Dirac operator V = 'y ■ {d + iA) + m: 

75P75 = (1) 

which, in particular, implies positivity of the determinant, detV > 0. For the correlator of 
a generic meson M = ip'^ip^ we can write, using (|l]) and Bunyakovsky-Schwartz inequality: 

{M{x)M^{0))^,A = -{TtS{x,0)TS{0,x)T)a 

= {TtS{x, 0)rt^,S\x, 0)ij5T)a < {TtS{x, 0)S\x, 0))a, (2) 

where S = T>~^ and V = 7or^7o. The inequality is saturated for mesons with T = i'j^Ti, 
since T> commutes with isospin r^, which means that the pseudoscalar correlators majorate 
all other 7 = 1 meson correlators.^] As a consequence, one obtains an important restriction 
on the pattern of the spontaneous symmetry breaking. For example, the symmetry breaking 
cannot be driven by a condensate of {ipj^ip). Indeed, broken axial SU(2) symmetry generators 
acting on such a pseudoscalar condensate would have produced O"*" Goldstone bosons ipTiip. 

At finite isospin density, fij 7^ 0, positivity still holds |]T5| and certain inequalities can 
be derived (in contrast with the case of /x^ 7^ when there is no positivity and hence no 
inequality can be derived). Now V = j ■ {d + iA) + |/i/7oT3 -|- m, and Eq. ([1|) is not true 
anymore, since the operation on the r.h.s. of (|l]) changes the relative sign of /i/. However, 

^It is important, as is the case for 1 = 1, that there is no disconnected piece after tp integration in (|^). 
The proof does not apply to cr- meson correlator, F = 1. 
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provided m„ = m^, interchanging up and down quarks compensates for this sign change (the 
u and d quarks play the role of mutually conjugate quarks fill), i.e, 



n75^^75ri = VK (3) 

Instead of isospin ti in one can also use T2 (but not r-^.) Eq. (H) replaces the now invalid 
Eq. (0) and ensures that detD > 0. Repeating the derivation of the QCD inequalities using 
(0) we find that the lightest meson, or the condensate, must be in channels 4'ij5Ti^24', i.e., a 
linear combination of tt~ ~ u'y^d and 7r+ ~ d'y^u states. Indeed, as shown below, in both two 
analytically tractable regimes of small and large /ii the lightest mode is a massless Goldstone 
mode which is a linear combination of ^75^ and d'j^u. 



3 Small isospin densities: pion condensate 

When hj is small, chiral perturbation theory can be used to treat the problem. To have 
a rough sense of how small /i/ should be, we require that no particles other than pions 
are excited due to the chemical potential. This gives /i/ ~ rrip as the upper limit of the 
applicability of the chiral perturbation theory. 

For zero quark mass and zero the pions are the massless Goldstone bosons of the 
spontaneously broken SU(2)iXSU(2)/j chiral symmetry. In reality, the quarks have small 
masses, which break this symmetry explicitly. Assuming equal quark masses, the symmetry 
of the Lagrangian is SU(2)2,+h- The low-energy dynamics is governed by the familiar chiral 
Lagrangian, which is written in terms of the matrix pion field S G SU(2): 

£ = ^f^TT[d,Ed^J:^ - 2m^ReE]. 

This Lagrangian contains only two phenomeno logical parameters: the pion decay constant, 
and the pion mass in vacuum, m-„. We will see that interesting physics occurs at /i/ > m^, 
and since m^r -C mp, there is a nontrivial range of /i/ where the chiral Lagrangian is a reliable 
and useful treatment. 

The isospin chemical potential further breaks SU(2)l+/j down to \J{1) l+r- Its effect 
can be included into the effective Lagrangian to leading order in /i/, without introducing 
additional phenomenological parameters. Indeed, /i/ enters the QCD Lagrangian in the 
same way as the the zeroth component of a gauge potential [|13|. Thus the finite- /i/ chiral 
Lagrangian is obtained by promoting the global SU(2)j;^x 811(2)/;; symmetry to a local gauge 



symmetry: gauge invariance completely fixes the way fij enters the chiral Lagrangian [13]: 

C2 ^2 f2 



C,s = ^TrV.SV.St - ^:^ReTrS. (4) 

The covariant derivative is defined as 

VoS = 9oS-^(r3S-Er3), V,S = (5) 

which follows from the transformation property of S under rotations by the isospin generator 
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Using (|^) it is straightforward to determine vacuum alignment of S as a function of /i/ 
and the spectrum of excitations around the vacuum. We will be interested in negative /i/, 
which favors neutrons over protons, as in neutron stars. The results are very similar to the 
two-color QCD at finite baryon density |jl3[- From (0), one finds the potential energy for S, 



V;ff(S) = :^Tr(r3Sr3St - 1) - ^^ReTrE . (6) 
o z 

The first term in (^) favors directions of S which anticommute with t^, i.e., ti and T2, while 
the second term prefers the vacuum direction E = 1. It turns out that the minima of (|^) at 
all Hi are captured by the following Ansatz: 

E = cos a + i{Ti cos (f) + T2 sin cj)) sin a . (7) 

Substituting to (^, one sees that the potential energy depends only on a, but not 0: 

Vesia) = ^(cos2a - 1) - /Xcosa . (8) 

Minimizing V^ff(a) with respect to a, one sees that the behavior of the system is different in 
two distinct regimes: 

(i) For l/i/l < m^, the system is in the same ground state as at /x/ = 0: a = 0, or E = 1. 
This result is easy to understand. The lowest lying pion state costs a positive energy 
— 1/^7 1 to excite, thus at zero temperature no pion is excited. The ground state of the 

Hamiltonian at such hj coincides with the normal vacuum of QCD. The isospin density is 
zero in this case. 

(ii) When \hi\ exceeds m.^ the minimum of occurs at 

2 

COS a = — (9) 

In this regime the energy to excite a tt" quantum, 7717^ — l/x/l, is negative, thus it is energetically 
favorable to excite a large number of these quanta. Since pions are bosons, the result is a 
Bose condensate of vr". If the pions did not interact, the density of the condensate would 
be infinite. However, the repulsion between pions stabilizes the system at a finite value of 
the isospin density. This value can be found by differentiating the ground state energy with 
respect to /xj: 

n: = = L^^JS.n a = f^^^J ^1 - ^ j . (10) 

For l/i/l just above the condensation threshold, — m^r <^ m,r, Eq. ([T0|) reproduces the 



equation of state of the dilute non-relativistic pion gas |l3 

= 4/^(/i/ - m^) 
At larger fij, ^ m^r, the isospin density is linear in /xj, 

nj = flixi, l/x/l > 
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From Eq. (10) one can find the pressure and the energy density as functions of /i/. The 
interesting quantity is the ratio between the two, 



e fxj + 3ml 



(11) 



This ratio starts from at threshold and quickly approaches 1 as one increases hj. Thus, as 
far as the chiral Lagrangian is still applicable, the EOS approaches the Zel'dovich limit of 
maximal stiffness at high densities, similar to nuclear matter in Walecka model. 

The fact that the minimum of the potential is degenerate with respect to the angle 
(j) corresponds to the spontaneous breaking of the 11(1)1,+^^ symmetry generated by in 
the Lagrangian (^). This is not unexpected since the ground state is, in essense, a pion 
superfluid, with one massless Goldstone mode. Since we start from a theory with three 
pions in vacuum, in the superfluid there are, in addition to the massless mode, two massive 
modes. One can be identified with 7r°. The other is a linear combination of vr"*" and vr", 
which we denote as tt~^, since it coincides with tt"*" at the condensation threshold. The mass 
(defined as the rest energy) of these modes can be obtained by expanding the Lagrangian 
(|) around the minimum. The result reads (cf. [p!3| ) 



At the condensation threshold, m^o = and m^f+ = 2mT,, while for m^r both masses 

approach (see Fig.|l|). 

The values of the chiral condensate, {uu + dd), and the pion condensate, (uj^d) follow 



i.e., the chiral condensate "rotates" into the pion condensate as a function of 

It is also possible to find baryon masses, i.e. the energy cost of introducing a single baryon 
into the system. The most interesting baryons are those with lowest energy and highest 
isospin, i.e. neutron n and A~ isobar. There are two effects of fij on the baryon masses. The 
first comes from the isospin of the baryons, which effectively reduces the neutron mass by 
and the A~ mass by If this was the only effect, the effective A~ mass would 

vanish at = |mA. For larger fij, baryon or antibaryon Fermi surfaces would form, which 
lead to a nonzero baryon susceptibility xb = dns/dfiB- However, long before that another 
effect turns on: the 7r~'s in the condensate tend to repel the baryons, lifting up their masses. 

These effects can be treated in the framework of the baryon chiral perturbation theory 
p!6| . For example, the (Euclidean) Lagrangian describing nucleons and their interactions 
with the pions at finite fij can be written as: 




(12) 



from d^): 



{uu + dd) = 2 ('?/;'?/') vac cos a and (u'j^d) + h.c. = 2 (■?/'■?/') vac sin a. 



(13) 




(14) 



where 
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Figure 1: Schematic plot of masses (rest energies) of lowest lying excitations in QCD at finite 
(negative) /ii, in the regime of applicability of chiral perturbation theory: ttIt^, /ij <^ nip. 
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Diagonalizing this bilinear Lagrangian in the pion background given by S = E from (J^) one 
finds the nucleon masses. The result for the neutron and the A~ isobar reads: 

l/^/l 3|/i/| 
m„ = mjv —cos a, m^- = m/^ — cos a, (15) 



in the approximation of nonrelativistic baryons. Equation (|T5|) can be interpreted as follows: 
as a result of the rotation (|^) of the chiral condensate, the nucleon mass eigenstate becomes 
a superposition of vacuum n and p states. The expectation value of the isospin in this state 
is proportional to cos a appearing in (15). With cos a given in Eq.(^, we see that the two 
mentioned effects cancel each other when m,r |yU/| -C rrip. Thus the baryon mass never 
drops to zero, and = at zero temperature in the region of applicability of the chiral 
Lagrangian. 

As one forces more pions into the condensate, the pions are packed closer and their 
interaction becomes stronger. When fij ~ rUp, the chiral perturbation theory breaks down. 
To find the equation of state in this regime, full QCD has to be employed. As we have 
seen, this can be done using present lattice techniques since the fermion sign problem is not 
present at finite /i/, similar to the two-color QCD p^ . 



4 Asymptotically high isospin densities: quark-anti- 
quark condensate 

In the opposite limit of very large isospin densities, or \fij\ ^ nip, the description in terms 
of quark degrees of freedom applies since the latter are weakly interacting due to asymptotic 
freedom. In our case of large negative /i/, or nj, the ground state contains an equal number 
of d quarks and u antiquarks per unit volume. If one neglects the interaction, the quarks fill 
two Fermi spheres with equal radii |/i/|/2. Turning on the interaction between the fermions 
leads to the instability with respect to the formation and condensation of Cooper pairs, 
similar to BCS instability in metals or the diquark pairing at high baryon density 0. To 
the leading order of perturbation theory, quarks interact via the one-gluon exchange. It is 
easy to see that the attraction is strongest in the color singlet channel, thus the Cooper pair 
consists of a u and a d. The ground state, hence, is a fermionic superfiuid. 

The perturbative one-gluon exchange, however, does not discriminate between the scalar, 
ud, channel, and the pseudoscalar, u'j^d, channel: the attraction is the same in both cases. 
But one can expect that the instanton-induced interaction, however small, will favor the 
u'-f^d channel over the ud one. The condensate thus is pseudoscalar and breaks parity, 

(u-f^d) ^ 0. (16) 

This is consistent with our earlier observation that QCD inequalities constrain the 1 = 1 
condensate to be a pseudoscalar at any /i/. Note that the order parameter in (^B[) has the 
same quantum numbers as the pion condensate at lower densities. We shall discuss this 
coincidence later. 

As a consequence of the Cooper pairing, the fermion spectrum acquires a gap A at the 
Fermi surface, where 

A = c = 37rV2, (17) 
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where g should be evaluated at the scale \fij\. The peculiar e~'^^^ behavior comes from 
the long-range magnetic interaction, as in the superconducting gap at large fiB [0]- The 
constant c is smaller by a factor of a/2 compared to the latter case due to the stronger one- 
gluon attraction in the singlet qq channel compared to the 3 diquark channel. Consequently, 
the gap (0) is exponentially larger than the diquark gap at comparable baryon chemical 
potentials. Using the methods of one can estimate b ~ 10^. As in the BCS theory, the 
critical temperature, at which the superfluid state is destroyed, is of order A. 

Asymptotically, A is much less than fij, and superfluidity has little effect on the equation 
of state. The ratio p/e approaches 1/3 from below in the limit fij — >■ oo. 



5 Quark-hadron continuity and confinement 



Since the order parameter ([T6|) has the same quantum numbers and breaks the same sym- 
metry as the pion condensate in the low-density regime, it is plausible that there is no phase 
transition along the fij axis. In this case, as one increases the density, the Bose condensate 
of weakly interacting pions smoothly transforms into the superfluid state of ud Cooper pairs. 
The situation is very similar to that of strongly-coupled superconductors with a "pseudogap" 
P^ , and possibly of high-temperature superconductors |]TB|. This also parallels the conti- 
nuity between nuclear and quark matter in three-flavor QCD as conjectured by Schafer and 
Wilczek We hence conjecture that in two- flavor QCD one can move continuously from 



the hadron phase to the quark phase without encountering a phase transition. We stress 
here that this conjecture needs to be verified by lattice calculations. 

At first sight, this conjecture seems to contradict a common wisdom that there is a 
"deconfinement" phase transition from the hadron phase to the quark phase. It is logically 
possible that there exists a first order phase transition at intermediate value of /i/. However, 
there are several nontrivial arguments that make the continuity hypothesis highly plausible. 

The first argument arises from considering baryons. One notices that all fermions have 
a gap at large which means that all excitations carrying baryon number are massive. 
In particular, at zero temperature, the baryon number susceptibility xb vanishes. This is 
also true at small fij. It is thus natural to expect that all excitations with nonzero baryon 
number are massive at any value of fij, and xb remains zero at T = for all fij. This also 
suggests one way to verify the continuity on the lattice. 

Another argument comes from considering the limit of large number of colors Nf.. Let 
us recall that in finite-temperature QCD, there is a mismatch, at large Nc, between the 
number of gluon degrees of freedom, which is 0{N^), and of hadrons, which is 0{N^). This 
fact is a strong hint of a first order confinement- deconfinement phase transition, at which 
the effective number of degrees of freedom jumps from 0{N^) to 0{N^). It is easy to see, 
however, that the A^^^ behavior of thermodynamic quantities is the same in the "hadronic" 
phase (low /i/) and the "quark" phase (large /i/). Indeed, at very large fij the isospin density 
rii is proportional to the number of quarks, which is 0{Nc): 

nj - 



3 87r2 

In the small /i/ region the isospin density is given by Eq. (plQl). At large Nc limit the pion 
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decay constant scales as = 0{Nc), and thus the isospin density in the pion gas is also 
proportional to N^.^ What happens is that the repulsion between pions becomes weaker as 
one goes to larger Nc, thus more pions can be stacked at a given chemical potential. As a 
result, the Nc dependence of thermodynamic quantities is the same in the quark and the 
hadronic regimes, although for seemingly very different reasons. 

Now let us return to the question of confinement. Naively, one would think that at 
asymptotically large fij, the u and d quarks are packed at a very high density, and the system 
should become deconfined. At finite temperature, there is no rigorous way to distinguish 
the confined and deconfined phases in QCD with quarks in the fundamental representation. 
However, at zero temperature (and finite /i/), a sharp distinction can be made between the 
two phases. In the confined phase, all particle excitations carry integer baryon number; 
the deconfined phase can be defined as the phase where there exist finite-energy excitations 
carrying fractional baryon charge. The pion superfluid at small /// clearly is in the confined 
phase. The question is: is the quark matter at large /i/ confined or deconfined? 

It might seem that at very large fij there exist excitations with fractional baryon number. 
These are the fermionic quasiparticles near the Fermi surface, which are related to the original 
quarks and antiquarks by a Bogoliubov-Valatin transformation. The opening of a BCS gap 
makes the energy of these excitations larger than A, but still finite. 

To see that the logic above has a fault and there are no such excitations, one needs to 
consider dynamics of very soft gluons. The crucial observation is that at large /i/, gluons 
softer than A are not screened neither by Meissner nor by Debye effect. Q Meissner effect 
is absent because the condensate does not break gauge symmetry (in contrast to the color 
superconducting condensate [Q). Debye screening is also absent, because on scales softer 
than A there are no charge excitations in the medium: the Cooper pairs are neutral, while 
the fermions are too heavy to be excited. Thus, the gluon sector below the A scale is 
described by pure SU(3) gluodynamics, which is a confining theory. This means there are 
no quark excitations above the ground state: all particles and holes must be confined in 
colorless objects, mesons and baryons, just like in vacuum QCD. 

If there is no transition along the fij axis, we expect confinement at all values of fij. At 
large fii, since the running strong coupling as at the scale of A is small, the confinement 
scale Aq(jj3 is much less than A. In more detail, let us imagine following the running of the 
strong coupling from the UV to the IR. First, as increases until the scale g^i is reached when 
it "freezes" due to Debye screening and Landau damping. The freezing continues until we 
reach the scale A, after which the coupling runs again as in pure gluodynamics. Since the 
coupling is still small at the scale A, it can become large only at some scale Aqqj-, much lower 
than A. Thus, at large \fii\ there are three different scales separated by large exponential 
factors, /ij > A > ^'qcd- 

That the scale of confinement is much smaller than the gap at large fij has an important 
consequence for finite temperature. One can actually predict a temperature driven deconfine- 

^ With physical values of Nc, fir and m7r, the values of ni given by eqs. (p^ ) and naively continued 
into the regime of intermediate /i/, cross at /x/ « SOOMeV. This agrees with the value of ^ rup where one 
would expect the crossover between the quark and hadron regimes to occur. This is a quantitative indication 
that a phase transition is not necessary. 

This is similar to the behavior of the unbroken SU(2)c sector of two-flavor color superconductors [ pO| . 
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ment phase transition at a temperature of order Aq^j-,. Indeed, at such low temperatures, 
quarks are unimportant, so the transition must be of first order as in pure gluodynamics. In 
particular, one expects the baryon number susceptibility temperature dependence to change 
from e~^^/-^ to e~^/^ around due to deconfinement. 

The smallness of the confinement scale Aqqj-, compared to the BCS gap A allows one 
to conclude that the binding energy of quarks and antiquarks is small and the hadronic 
spectrum follows the pattern of the constituent quark model, with A playing the role of the 
constituent quark mass. This means mesons weigh 2A and baryons weigh 3A, approximately. 
A good analog of the large iij regime is vacuum QCD with only heavy quarks. As in the 
latter case, the string tension and string breaking are determined by parametrically different 
energy scales (Aqq^ and A, respectively). Hence the area law should work up to some 
distance much larger Aq^q, even when fundamental quarks are present. For the same reason 
one also expects the high-spin excited states of hadrons to be narrow at large fij. 

The energy hierarchy also leads to a curious dispersions relation of hadrons in the isospin 
dense regime. Consider, for example, the p~ meson, which is a bound state of a m and a d. 
At zero total momentum, the u and the d are on the opposite sides of the Fermi surface. 
As one increases the total momentum, the two constituents move along the Fermi surface, 
remaining close to the latter until the total momentum becomes larger than /^/, i.e. twice 
the Fermi momentum. Thus, the dispersion curve of the p~ must remain essentially flat 
in the interval of momentum (0, For baryons, energy is almost independent of the 

momentum in the interval (0, 1.5|/i/|). The group velocity of hadrons, thus, almost vanishes 
in these intervals. Above these intervals it should be equal to the speed of light. It would be 
interesting to follow, on the lattice, the evolution of the dispersion curves of p~ from small 
to large pj. 

6 The phase diagrams on (T, fij) and (/i, /ij) planes 

By considering nonzero pi, we make the phase diagram of QCD three-dimensional: (T, ps, 
Pi). Two planes in this three-dimensional space are of a special interest: the Pb = (T,pj) 
plane, which is completely accessible by present lattice techniques, and the T = {pi.pb) 
plane, where the neutron star matter belongs. 

Let us first consider the phase diagram on the (T, pi) plane, which is simpler. Two phe- 
nomena determine the phase plane on this plane (Fig.^: pion condensation and confinement. 
At sufficiently high temperature the condensate (|16D melts (solid line in Fig. |^). For large pi, 
this critical temperature is proportional to the BCS gap (|T7p. There are two phases which 
differ by symmetry: the high temperature phase where the explicit flavor U{1)l+r symmetry 
is restored, and the low-temperature phase where this symmetry is spontaneously broken. 
The phase transition is in the 0(2) universality class.0 The critical temperature vanishes 
at Pi = and is an increasing function of pi in both regimes we studied: \pi\ <^ nip and 
\pi\ ^ Aqcd- Thus, it is likely that Tc{pi) is a monotonous function of pi. In addition, as 
explained before, at large /i/, there is a first-order deconfinement phase transition at some 

* The width of the Ginzburg region is suppressed by (A//i/)^ at large /i/ as in usual BCS superconductors 
and also by l/iV^ at large Nc as at the QCD chiral transition [ pT] . 
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Figure 2: Phase diagram of QCD at finite isospin density. 



temperature much lower than Tc{fii). Since there is no phase transition at /i/ = (for 
small m^rf) or at T = (assuming quark-hadron continuity), this first-order line must end 
at some point A on the {T^jij) plane (Fig. 

The phase diagram in the (/x/, /i) plane at zero temperature turns out to be quite com- 
plicated. We defer the more detail study of this plane for future work. Here we shall only 
consider the regime |yU/| ^ fiB, both much larger than Aqcd, so that perturbative QCD 
can be used. When /i^ = and ^ Aqcd we have seen that the system is a superfluid 
with a gap A. Finite fiB provides a mismatch between u and d Fermi spheres. The su- 



perconducting state becomes unfavorable at some value of /ig of order A. It is known p2 



that the destruction of this state occurs through two separate phase transitions. As one in- 
creases hb, at hb slightly below A/ \/2, a first-order phase transition takes the system to the 
Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) state |2^, which is characterized by a spatially 
modulated superfluid order parameter (u'jQd) with a wavenumber of order 2fiB- How exactly 
this spatial dependence looks like is still unknown, mostly because the FFLO state has not 
been observed in metals. The FFLO state persists only until fiB = 0.754A when it goes 
through a second-order phase transition to a {u'j^d) = state. The latter must be a color 
superconductor with one- flavor diquark condensates (uu) and {dd), due to the attraction 
between quarks of same flavor. In the region of the (fiBfJ^i) phase diagram directly relevant 
to neutron stars, fiB > fii, the color superconducting FFLO phase is studied in a recent 
paper of Alford, Bower and Rajagopal 



The most interesting feature of the FFLO state is that it has the same symmetries as 
the inhomogeneous pion condensation state which might form in electrically neutral nuclear 
matter at high densities, as argued by Migdal and others 0. The FFLO phase, thus, can 
be thought of as a realization of Migdal's pion condensate in the regime of asymptotically 
high densities. It is also conceivable that the two phases are actually one, i.e., continuously 
connected on the {fij,fiB) phase diagram. 



12 



7 Conclusion 



Our original and primary motivation for considering QCD at finite isospin densities is to 
fiave a dense regime of realistic, three-color QCD that can be studied on lattice. Based on 
analytical calculations in the two asymptotic regimes of low and high densities, we found 
that there is likely no phase transition along the /i/ axis at zero temperature. This conjecture 
should be verified on the lattice. An obvious way is to study the thermodynamics of the 
system. If our continuity conjecture is correct, all thermodynamic quantities should be 
smooth functions of fij. In this case, we also suggest that the ratio p/e is a non-monotonic 
function of /i/: it raises from at the threshold /i/ = 171^^ to some value close to 1, then 
drops to some minimal value, and then approaches 1/3 from below at large fij. The baryon 
susceptibility should vanish at any fij at zero temperature. We also predict the existence of 
a line of first-order phase transition on the (T, /i/) plane, which terminates at a second order 
point. 

The phase diagram on the /xb) plane, which is most relevant for neutron star physics, 
remains inaccessible to the lattice. Based on our preliminary investigations, the phase dia- 
gram on this plane should have a rather complicated topology. The most interesting feature 
of this diagram appears to be the existence of the FFLO phase, which is reliably predicted 
at iJ,i ^ fiB, both being large. This phase has the same symmetry as the pion condensation 
state conjectured by Migdal, and both might be different regions of a single connected region 
on the phase diagram. 

The authors thank L. McLerran, J. Kogut, R. Pisarski, E. Shuryak for discussions, K. Ra- 
jagopal for drawing their attention to Ref. ||2^, and the DOE Institute for Nuclear Theory 
at the University of Washington for its hospitality. The work of DTS is supported, in part, 
by DOE grant DE/FG02-92ER40699. 



Note Added 



Further research reveals that the physics below the scale A is described by "gluodynamics 
of continuous media" with a large dielectric constant. As a result, the scale of the confinement 
^QCD is small and decreases exponentially with the chemical potential fij. This means that 
the line of the first order deconfinement transition goes down as indicated in Fig. & 
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